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Motivated by recent experiments on SrQi2(B03)2, we investigate the ground states of the doped Mott insula- 
tor on the Shastry-Sutherland lattice. To provide a unified theoretical framework for both the valence-bond solid 
state found in undoped SrCu2(BC>3)2 and the doped counterpart being pursued in on-going experiments, we an- 
alyze the t-J-V model via the bond operator formulation. It is found that novel superconducting states emerge 
upon doping with their properties crucially depending on the underlying valence bond order. Implications to 
future experiments are discussed. 
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I. INTRODUCTION 

The fate of doped Mott insulators is one of the most chal- 
lenging issues in correlated electron physics, especially in re- 
lation to the long-standing problem of high-temperature su- 
perconductivity in cupratesii 2 ^ 3 -. Recent discoveries of vari- 
ous Mott insulators on geometrically frustrated lattice o 4 ' 5 ' 6 ' 7 
and organic materials 8 may offer an important clue to this is- 
sue when such materials are doped. While it is possible that 
the ground state does not break any symmetry, resulting in 
a spin liquid phase, the ground states of Mott insulators of- 
ten have broken spin-rotation and lattice-translation symme- 
tries, leading to antiferromagnetic and valence bond solid or- 
der, respectively-. Since both antiferromagnetic and valence 
bond solid phases are generic possibilities for Mott insulators, 
a zero-temperature quantum phase transition may occur be- 
tween these two phases when an appropriate "control param- 
eter" is changed 1 . In this context, understanding the effect of 
doping on valence-bond solid phases is as equally important 
as that on antiferromagnetic phases and would be quite useful 
for the full classification of all possible phases of doped Mott 
insulators. 

There are, however, not many clear examples of two- 
dimensional valence-bond solid insulator in contrast to the 
antiferromagnetic insulator found in high T c cuprate com- 
pounds. The discovery of SrCu 2 (B03)2 is particularly im- 
portant in this regar d 4 ' 9 ' 10 . This material can be charac- 
terized by an antiferromagnetic spin- 1/2 Heisenberg model 
on the Shastry-Sutherland lattice^. Starting from the usual 
square lattice, the Shastry-Sutherland lattice can be obtained 
by putting additional diagonal bonds with two possible ori- 
entations in alternating plaquettes (see Fig.l). Let J and J' 
be the exchange couplings along the diagonal and the square 
lattice links, respectively. It is known that the valence bond 
solid state or the product state of valence bond singlets (illus- 
trated as filled ellipses in Fig.l) along the diagonal bonds is 
the exact ground state of the Heisenberg model for J 1 / J < 
0.7 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 . In SrCu 2 (B0 3 ) 2 , J' /J is estimated to 
be 0.641. 

In this paper, we investigate possible phases of doped 
valence-bond solid and superconductivity on the Shastry- 
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FIG. 1: (Color online) Schematic diagram for the ground states on 
the Shastry-Sutherland lattice. The valence-bond solid state is the 
ground state at half filling as depicted in the left figure where ellipses 
represent spin-singlet pairs. When doped with holes, the system ex- 
hibits superconductivity in addition to the coexisting valence-bond 
solid order. Furthermore, when the nearest-neighbor repulsive inter- 
action, V, is larger than a critical value, V c , the plaquette D-wave su- 
perconductivity appears in a range of doping concentration, x, with a 
peculiar spatial pattern as shown in the right figure. In this situation, 
the pairing amplitudes residing in the four links encircling the hori- 
zontal dimers have the opposite sign to those for the vertical dimers. 
Different colors are used to emphasize the plaquette pattern of the 
pairing amplitude. 



Sutherland lattice. In order to study the interplay between 
the valence-bond solid order and emergent superconductivity 
at finite doping, we use the bond operator formulatio n 19 ' 20 ' 21 
of the constrained Hilbert space of correlated electron sys- 
tems, extended to general doping and applied to the t-J-V 
model. Here t and V represent the hopping strength and 
nearest-neighbor repulsion between electrons, respectively. 
In contrast to previous studies of the Shastry-Sutherland 
model 23,24 ' 25 , the emergent superconducting state is directly 
related to the underlying valence-bond solid order at the half- 
filling. For example, the valence-bond solid order is so robust 
that superconductivity always coexists with it as shown in the 
resulting phase diagram of Fig. 2. 

When the nearest-neighbor repulsion, V, is smaller than a 
critical value V c , the doped holes can remian paired mostly 
within each valence-bond singlet or dimer on the diagonal 



2 



1.5 
V/t ! 

0.5 



\laquet e 
i-wav< 



-Mott-insulatori 



S-wave 



S+\D wave - 



Normal Metal 



-Vc 



0.25 0.5 0.75 1 

X 

FIG. 2: (Color online) Phase diagram of the t-J-V model on the 
Shastry-Sutherland lattice as a function of hole concentration, x, and 
the nearest-neighbor repulsive interaction, V. The thick line sepa- 
rating the normal metal and the S-wave superconducting phase is 
a first-order phase boundary, while other boundaries are all second- 
order. The undoped Mott-insulating state is represented by the thick 
solid line at x = 0. 



bonds. This leads to S'-wave superconductivity. On the other 
hand, when V is larger than V c , the Coulomb repulsion pre- 
vents doped holes from occupying the same dimer. At larger 
doping concentrations, the doped holes tend to develop pair- 
ing correlation between nearby dimers of the valence-bond 
singlets. Consequently the superconducting order paramter 
acquires an interesting phase relation on the square lattice 
link, depending on whether the square lattice links encircle 
the horizontal or vertical dimers (See Fig. 1). More specif- 
ically, all the links within the same plaquette have the same 
sign, but those in the nearest-neighbor plaquettes have the op- 
posite sign. Since this looks like D-wave symmetry at long 
distances, we call it the plaquette D-wave superconductiv- 
ity. Notice, however, that this is different from the ordinary 
d-wave superconductor on the square lattice where the sign 
alternates for different links within the same plaquette. This 
peculiar structure comes from the underlying valence bond 
solid order. The simultaneous presence of the valence bond 
order and superconductivity can, in principle, be checked by 
X-ray or neutron scattering experiments. 

The rest of the paper is organized as follows. In section 
HH we review the bond operator formulation for doped mag- 
nets. In section [TTTJ the t-J-V Hamiltonian is written in the 
bond operator representation and the mean field theory is ex- 
plained. In section |IV| the results of the mean field theory 
and the phase diagram are discussed. We conclude in section 
[Vl Some details of the computations are relegated to the ap- 
pendix. 



n. BOND OPERATOR FORMALISM 

We start with the bond operator theory by setting up an ex- 
act mapping between bond operators and usual electron cre- 
ation operators 19 i22i2i. Let c\ a and c\ a (a =f, J.) be the elec- 
tron creation operators on the left and right site, respectively, 
of the lattice link that a horizontal dimer occupies (see Fig.Q]). 



In the limit of large on-site repulsive interaction, U, any state 
with two electrons at the same site is excluded from the low- 
energy Hilbert space which in turn can be represented by nine 
"bond-particle" creation operators defined as follows: 
(i) Singlet boson for spin-Peierls order, 

s^\v) = -ie«6cL<4|0>, (1) 



where |0) is the electron vacuum and \v) is the imaginary vac- 
uum void of any bond particle. E a j, is the second-rank anti- 
symmetric tensor with e-\i = 1. From now on we adopt the 
summation convention for repeated spin indices (such as a and 
b above) throughout this paper unless mentioned otherwise. 

(ii) Triplet magnon, 

4» = ^<EcacL4|0), (2) 

where a" b (a — x, y, z) are the usual Pauli matrices. 

(iii) Single-hole fermion, 

4» = 4lo>, 

4»=4a|0>- (3) 

(iv) Double-hole boson for the empty state, 

S\v) = |0). (4) 

As indicated by the names, the operators s, d, and t a all 
obey the canonical boson commutation relations while h\ a 
and 1%2a satisfy the canonical fermion anticommutation rela- 
tions. To eliminate unphysical states from the enlarged Hilbert 
space, the following constraint needs to be imposed on the 
bond-particle Hilbert space: 



s's ■ 



tU 



h\ a hia + h\ a h 2a + d)d = 1. 



(5) 



Constrained by this equation, the exact expressions for the 
spin and electron creation operators can be written in terms 
of the bond operators. For example, 



S 2 a 



Aa = h la d + ^£afcS f /l2h - ^=E ac a%t\h 2b , 



J 
'2 a 



hLd + 
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SacVcbAxhlb, (6) 



where e Q /3 7 is the third rank antisymmetric tensor with e xyz = 
1. 

Now, in order to cover the full Shastry-Sutherland lat- 
tice, we need nine additional bond operators representing the 
electronic Hilbert space associated with the vertical dimers 
corresponding to the sites 3 and 4. For example, r Q indi- 
cates the triplet magnon in the vertical dimers while h^ a and 
hi a represent the corresponding single-hole fermions. While 
non-uniform condensations of the spin-Peierls singlet and/or 
double-hole bosons are in general possible, we focus only on 
the uniform phases in this study. Thus the condensation densi- 
ties of the singlet and double-hole bosons will be represented 
by s and d for both horizontal and vertical dimers. 
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III. CHOICE OF HAMILTONIAN AND BOND OPERATOR 
MEAN FIELD THEORY 



We consider the following t-J-V Hamiltonian: 



H = Pg[-J2 *« ^ia C ia + H. C.) ~ M 2 C i* Cia 

+ ^2 ^i Sj ' S i + v v n i n 3~\ Pg , (7) 

(hi) (*.i> 



where Pc is the Gutzwiller projection operator imposing the 
no-double-occupancy constraint. Uj and Jy are the electron 
hopping matrix element and antiferromagnetic exchange in- 
teraction, respectively. In this paper, we set £y = t and 
= J within a dimer, tij = t' and = J' between neigh- 
boring dimers. We take J' / J — 0.64 from the experiments 
on SrCu2(B03)2 4,23 . From the large-C/ expansion of the Hub- 
bard model, we can then set t'/t = J J' / J = 0.8. However, 
the parameter, J'/t', which is important for hole dynamics, is 
not available experimentally. Thus we follow the convention 
used in previous works^ and take J'/t' = 0.3. Also, [i is 
the chemical potential and n, is the electron density operator. 
Finally, for convenience, the nearest-neighbor repulsive inter- 
action ,Vij, is set to V for both the two sites within dimer and 
those between the nearest-neighbor dimers. 

We now write the t-J-V Hamiltonian solely in terms of 
the bond particle operators 19 . As usual, the constraint on the 
bond particle operators is imposed by the Lagrange multiplier 
method. Residual interactions between bond particles are an- 
alyzed via quadratic decoupling of quartic terms in a simi- 
lar manner to the usual Hartree-Fock-BCS treatment. For the 
mean-field description, we consider the following order pa- 
rameters: 



Px — (j*i a 1~i-\-x,at) ; Q x — (tia1~i+x ,a) : 

Tlx = (h lia h3 } i-x,a) — (h>u a h4,i-x,a) 

= (^2ia^34+i,a) = (^iaVi+fi.a) i 

A x = (/llix/l3,i-x,T) = (^lil^4,i-S,T) 

= (^2ii^34+i,t) = (^2iJ,^4,i+*,t)' 



(8) 



where i is the dimer index of the horizontal dimers and i ± x 
indicates the locations of the neighboring vertical dimers. The 
order parameters for the y-direction can be defined similarly. 

We consider the possibility that the valence-bond-solid or- 
der persists even at non-zero doping. We further consider the 
condensation of d bosons, but neglect the possibility of triplet 
condensation because we are mostly concerned about param- 
agnetic phases in this work. All of the order parameters (P, Q, 
n, and A for both x and y directions) defined above as well as 
s, d, the chemical potential, fj,, and the Lagrange multiplier, £, 
are determined by solving a coupled set of four saddle-point 
equations and eight self-consistency equations. Details of the 
computational procedure are presented in the Appendix. 
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FIG. 3: (Color online) Superconducting order parameter, x> as we H 
as condensation densities of the spin-Peierls singlets, s, and double- 
hole bosons, d, as a function of hole concentration, x. When V = 
OAt > V c (top panel), S-wave \ an d d simultaneously vanish at x = 
0.431. Before the full plaquette D-wave superconductivity sets in at 
larger x, there is a region, 0.389 < x < 0.431, where the S- and D- 
wave superconductivity coexists in the form of S+iD-wave. Finally 
at x = 0.5, the plaquette D-wave superconducting order parameter 
becomes zero and so does s. No valence-bond solid correlation exists 
for x > 0.5. When V = O.lt < V c (bottom panel), d increases 
monotonically as a function of x while s remains finite all the way 
to a; = 1. For clarity the S- and plaquette D-wave superconducting 
order parameters as well as d 2 are magnified twenty times in the top 
panel while, in the bottom panel, only the S-wave superconducting 
order parameter is magnified twenty times. 



IV. RESULTS OF THE MEAN FIELD THEORY AND THE 
PHASE DIAGRAM 



At zero doping, the system is in a robust valence-bond solid 
phase. In our bond operator theory, the robustness of the 
ground state can be seen by the complete localization of triplet 
magnons: the t a and r Q dispersions are completely flat and 
high in energy. It can easily be shown that, at the quadratic 
order, the triplet magnons on the Shastry-Sutherland lattice 
are completely decoupled from the singlet contributions (No- 
tice that the coupling between the singlets and triplets at the 
quadratic order were the main driving force for the triplet dis- 
persion in the square-lattice cas©i- 9 ). This leads to the triplet 
Hamiltonian with no dispersive quadratic part. Thus, the only 
way to generate the triplet dispersion is through the saddle- 
point order parameters, P x , P y ,Q x and Q y . It turns out, how- 
ever, that all of the above order parameters are actually zero 
for any J'/ J at half filling in our mean-field theory so that the 
triplets are completely localized. 

Now let us consider the case of non-zero doping. In our 
Hartree-Fock-BCS saddle-point approximation, superconduc- 
tivity is found to appear at any non-zero doping when one 
neglects long-rage charge inhomogeneities such as Wigner 
crystal order at very small x. In the bond operator represen- 
tation, the electron superconducting order parameter can be 
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FIG. 4: /i-fermion pairing amplitude, A, as a function of hole con- 
centration, x. In the plot, As and Ad indicate the S- and D- 
wave components of A, respectively. Note that there is a region, 
0.389 < x < 0.431, where the S- and D-wave superconductivity 
coexist in the form of the S'+iD-wave. 



computed from 

Xx = e ab( c lia C 3,i-£,b) = £ ab ( C lia C 4,i-fi ,&) 

t J 



£ ab\C 2ia C 3i+£b ) — £ ab\C 2ia C 4i+ ~ b 



(2d 2 



Q X )A X - V2sdH x 



(9) 



The superconducting order parameter for the ^-direction is de- 
fined similarly. The doping dependence of \ is plotted in Fig. 
[3] along with the singlet boson condensate density, s, and the 
double-hole boson condensate density, d. 

While x can be a complex number in general, it turns out to 
be real for most of the phase space, barring an arbitrary overall 
phase factor (See Fig. [2]). In Fig. [2] however, there is a narrow 
region where the S- and D-wave superconductivity coexists 
in the form of an .S+iD-state similar to the results of previous 
studies on the square lattice 28 . In this case Xx is a complex 
number and Xx = X* y - As shown in Eq. (O, the /i-fermion 
pairing amplitude, A, is the key element in determining the 
sign of x an d thereby the pairing symmetry of the electronic 
superconducting order parameter. Generally the momentum 
dependence of A can be written as follows: 



As+iD = A x cos k x + A y cos k y 
= Ae iip cos k x + Ae 
= A s (cosk x - 



^ cos ky 
i Ad (cos 



- cos k y ) + 



k>£ COS ky) • 

(10) 



If both As and Ap remain finite, the system possesses the 
S'+iD-wave pairing symmetry. In Fig. [4] we plot the detailed 
doping dependence of As and Ajj for V = OAt. 

Notice that the critical exponent, /3, defined in x ~ \x— %c\ 
shows an interesting behavior when V > V c . Supercon- 
ducting order parameters for both S- and D-wave show the 
usual mean-field behavior, i.e., (3 — 1/2, when they first 
emerge. On the other hand, when they disappear, both ex- 
hibit an unusual exponent of j3 — 1. This deviation from the 
conventional mean-field behavior has different origins in each 
case. For the S-wave case, the phase boundary corresponds 



to the critical point between S+iD-wave and D-wave super- 
conductivity, where the presence of nodal fermions leads to a 
non-analytic cubic term in the expansion of the ground state 
energy 2 ^. On the other hand, for the D-wave, the linearly 
vanishing s 2 gives rise to j3 = 1 (Note that A d remains finite 
at the critical point x = 0.5. See Fig. [4]). 

It is also interesting to note that our /i-fermion (defined on 
a dimer) is an extended object which carries both the charge 
and spin quantum numbers. As a consequence, when the h- 
fermion pairing amplitude is finite, the distance between two 
holes in a h-h pair is larger than the average electron distance, 
as reported in exact diagonalization studies 29 . This behavior is 
different from what one would expect in the slave-boson-type 
theor y 23 ' 24 . 

In Fig. |2] the zero-temperature phase diagram is plotted as 
a function of hole concentration, x, and the nearest-neighbor 
repulsive interaction, V. The overall shape of the phase di- 
agram is determined by the behavior of the e?-boson conden- 
sate, which corresponds to the local pairing of holes within 
dimers. For V < V c ~ 0.23t, the d-boson condensation is the 
primary mechanism of pairing and its short-range nature gen- 
erates S'-wave pairing rather than D-wave with a nodal struc- 
ture. On the other hand, for V > V c , there is a V^-dependent 
critical hole concentration, x c , where d vanishes. It is this col- 
lapse of d-boson condensates that makes pairing more long- 
ranged and finally leads to the emergence of D-wave super- 
conductivity. 

The detailed ^-dependence of the superconducting order 
parameter, x, is shown in Fig. [3] For V = OAt > V c , the S'- 
wave superconducting order parameter and the rf-boson con- 
densate density go to zero simultaneously at x ~ 0.431. The 
collapse of the d-boson condensate at finite doping, in turn, 
imposes a precise upper bound on the hole concentration, up 
to which the valence-bond order can exist. To see this, note 
that, in the bond-operator representation, the chemical poten- 
tial, pt,, is fixed to satisfy (h\ ia hu a + h\ ia h-2\a) + 2d 2 = 2x. A 
similar equation exists for hs and fermions. When d 2 = 0, 
hole doping can be achieved only through /i-fermions. In 
this situation, the hole concentration cannot exceed x = 1/2 
where every dimer has exactly one electron. Note that, for 
V = OAt < V c , d 2 increases monotonically as x increases 
and finally reaches unity at x = 1. 

The plaquette D-wave superconducting order parameter 
found here has a peculiar spatial pattern as shown in Fig. Q] 
As explained in the introduction, this pattern is completely 
different from the corresponding pattern of the conventional 
d-wave superconducting state that previous slave-boson the- 
ories would predict 23,24 . We argue that the order parameter 
pattern of our D-wave state is rather natural in the limit of 
strong dimerization. In this limit, two electrons within the 
same dimer lose their independent local coordinates and share 
the same single dimer index. Then, the links of the plaquette 
enclosing the horizontal (vertical) dimers would connect only 
the links of the nearest-neighbor plaquettes that encircle the 
vertical (horizontal) dimers. Thus it is natural to expect the 
plaquette-dependent order parameter. 

Completely different behaviors of the d-boson condensate 
density below and above V c can physically be understood as 



5 



<€Z1> 




AE=E r E 2 



-2 

d 05 



V>V c (AE>0) 


1 : 




V<V c (AE<0)_ 



0.5 
X 



FIG. 5: (Color online) Schematic diagram showing the origin of dif- 
ferent behaviors of the d boson condensate density below and above 
V c ~ 0.23t. The bottom figure plots <P as a function of hole concen- 
tration, x, for V = OAt (AE > 0) and V = O.lt (AE < 0). Note 
that for clarity eP is magnified thirty times in the case of V = OAt. 



follows. When two holes are doped into a unit cell, there are 
two possible configurations as shown in Fig. [5] First, two 
electrons are removed from the same dimer while the other 
remains intact. In this case, only one singlet pair is removed, 
which costs the energy of E\ = E s ~ J with an additional 
Coulomb energy gain due to the disappearance of nearest- 
neighboring electrons: ~ V. Thus the total energy cost 
would be Ei = E s — Ed ~ J — V. Second, when two 
electrons are removed from two dimers, i.e., one hole remains 
at each dimer, there are twice the Coulomb energy gain and 
twice the exchange energy cost for breaking two singlet pairs. 
The total energy cost for the second case is E2 ~ 2E\. In 
conclusion, the energy-cost difference between two configu- 
rations is AE = E x - E 2 ~ V - J. If AE < 0, the first 
configuration is preferred and the d bosons can condense. On 
the other hand, if AE > 0, the second configuration is favored 
and the d bosons may not condense depending on x, which is 
precisely the case as shown in the bottom figure of Fig. [5] 



V. CONCLUSION 

By analyzing the t-J-V model based on the bond opera- 
tor formalism, we have investigated the phase diagram of the 
doped Mott insulator on the Shastry-Sutherland lattice. The 
interplay between strong dimerization and nearest-neighbor 
repulsive interactions leads to different behaviors of the doped 
holes determining the overall phase diagram. If the nearest- 
neighbor repulsive interaction, V, is smaller than a critical 
value, V c , the hole-pairing within dimers is preferred, result- 
ing in 5-wave superconductivity at any non-zero x. On the 
other hand, if V is larger than V c , the density of paired-holes 
within the same dimers vanishes at finite x and the plaque- 
tte Z?-wave superconductivity with a peculiar spatial pattern 
emerges. 

Considering clear experimental evidence for the valence- 



bond solid state in undoped SrCu2(B03)2, we believe that 
the above conclusions would be valid for a range of real- 
istic situations. It will be interesting to observe the possi- 
ble transition between the S- and D-wave superconducting 
states. Finally we note that, while doping mobile carriers to 
SrCu2(B03)2 has not been achieved yet, there has been a re- 
cent progress in doping quenched non-magnetic Mg impuri- 
ties to this material^. 

This work was supported by the NSERC, CRC, CIAR, 
KRF-2005-070-C00044 (YBK) and by the KOSEF through 
CSCMR SRC (BJY, JY). 



APPENDIX A: COMPUTATIONAL DETAILS 

Here we present details of the Hartree-Fock-BCS mean 
field theory in the bond operator formalism. We begin by ex- 
plicitly writing the t-J-V Hamiltonian as follows: 



H = P (H t + Hj + + H v ) P Q , 



(Al) 



where 



H t = -t^2(c\ ia c 2ia + c\ ia c Ua + H. c.) 

i 

^ ^ ^ |^lia (^3,i — x,a x,a) 
i 

+ 4ia( c 3,*+a,a + c 4.i+x.a) + H. C 

i 

-y,a + C2,i-y,a) + H. C. , (A2) 

Hj = J y^(Sij • S 2 ; + S 3i • S 4 j) 

i 

+ J' ' (^3,1-4 + ^i,i-x) + S2; ' (S3,i+^ + S^t+jfc 

i 

+ J' J]] [S4i • (Si,i+y + ^2,i+y) + &3i ' {&l,i-y + ^2,i-y 

(A3) 

Hp = -fi / ](nu + n 2l - 2 + 2.x) 

i 

- nJ2(n 3l + n il -2 + 2x), (A4) 

i 

and 

H v = V ^(niin 2l + n 3i n 4l ) 

i 

i 

+ V [ n 4*(ii,»-H/ + n 2 ,i+y) + n 3l (n M _y + n 2 ,i-y) 

(A5) 

In the next section, the t-J-V Hamiltonian is expressed solely 
in terms of bond particles. 
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1. Hamiltonian in the bond-operator representation 



and 



Taking usual steps in any saddle-point theory, we first re- 
place the full Gutzwiller projection by adding the Lagrange 
multiplier term. In the bond operator representation, this La- 
grange multiplier term is written as follows: 

Hi = -£^( S ! S * + 4a*i« + h\ ia hlia + h\ ia h 2% a + d\di - 1) 
i 

- £ y2( a i a i + T L r i<* + h La h 3*a + h\ la h 4l a + S}S t -l), 



Hi 



fen? 



(A6) 



where £ is the Lagrange multiplier. Also, er^ and Si represent 
the singlet boson and double-hole boson operators for the i- 
th vertical dimer, respectively. Other operators are defined in 
the main text. Under the constraint imposed by Eq. JA61 . the 
usual chemical potential term can be written as follows: 

H t* = ^J2( h Ua h Ha + h\ ia h 2la + 2djdi - 2 X ) 

i 

+ V X^ l 3i>3m + h\ la h Aia + 2djSi - 2x). (A7) 

i 

Since we focus only on the phases with homogeneous singlet 
boson and double-hole boson condensates, from now on we 
set s l = s\ = (7 [ = erf = s and d t = d\ = Si = Sj = d. 
Furthermore, we neglect the possibility of triplet condensation 
in this paper. 

By applying the bond operator representation to all the re- 
maining terms in the Hamiltonian, the saddle-point Hamilto- 
nian can be written as follows: 

H = Ne + fftriplet + -ff fermion, (A8) 

where the triplet boson Hamiltonian is given by 

-ff triplet = E ^k[*a (k)T a (k) + T+(k)t a (k)] 
k 

+ B k [tl (k)rl (-k) + r a (-k)t a (k)] 

k 

+ 5^C k [t+(k)t a (k) +rt(k)r a (k)], (A9) 

k 

and the single-hole fermion Hamiltonian is given by 



i rrinteraction / a t r\\ 

tl fermion — -"fermion fermion I \>^^) 



where 



H 



(0) 

fermion 



: J2 {a^h\ a (k)h 2a (k) + ht a (k)h 4a (k)\ + H. c.} 

k 

J2h[h{ a (k)h la (k)+hi(k)h 2a (k) 



hi(k)h 3a (k) + h\ a (k)h 4a (k) 



: E {c k [hl a (k)h la (k) + hl a (k)h 4a (k)} + H. c.} 

k 

E {d k [hl(k)h la (k) + hl(k)h 3a (k)} + H. c.} 

k 

E {eik[/4 T (k)fct ; (-k) + ^ T (k)ft^(-k)] +H. c.} 

k 

E { e 2k[^ 1T (k)^ ; (-k) + 4 T (k)4 ; (-k)] +H. c.} 

k 

E {M h h( k ) h n(-V + ^ T (k)^ x (-k)] + h. c.} 

k 

E {/2k[4 T (k)4;(-k) + 4 T (k)4 4 (-k)] +H. c.}. 



(A12) 



In the above expressions, N is the number of unit cells, 

3 



e = -^Js 2 + 2f (1 - t + d 2 - 2x) - 2 l i(2x - 2d 1 ) 



f + d 2 

J' (Ql + Ql~ Pi - Py) 



ij2 I OT/71 ~\2 



and 



+ V(2 + 2d 2 - Ax) + 8V(1 - x) 2 

A k = J'(P X cos k x + P y cos k y ), 

B k = -J'{Q X cosk x + Q y cos fey), 

C k = J/4-£, (A13) 



a k = -t, b k = fi, 

Ck = CRb + ici k , d k = d Rk + idm, 

eik = e Rk + ieik, e 2 k = eRk - ieik, 

hk = fnk + ifik, /2k = /flk - if ik, (A 14) 



where 



CRk = -t [d - -s ) (cosfer + COS fey) 

"— (P x cos k x + Py cos k y ) 
3 

- J'(H X cos k x + Uy cos fey), 



(All) 



cik = —t [d + -s ) (sinfc x - sinfer) 

"— (-P x sin fc x + P y sin fe y ) 
3 

- J' (He sin fc x - U y sin fey ) , 



(A15) 
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d Rk = -t' [d 2 



1 



s ) (cos fez + cos ky) 



(P x cos k 

x Py COS ky 



3 

- J'(Il x COS fca; + Ily COS A) y ), 



dj k = —t yd + —s ) (sin fca; + sin k v ) 



(P x sin fc^ + P y sin fc v 



3 

- J'(H X sink x + Il y sinky), (A16) 



e_Rk = — v2Y ds(cos fc x + cosk y ) 
3 

— — J' (A x cos fcj; + A a cos k y ), 
3 

e/k = -- J'(A X sink x - A y sinfc a ), (A17) 



and 



fm = — \/2t' ds(cos k x + cos k y ) 
3 

— -J'(A X cos k x + A y cos fcj,), 
3 

fik = --J\A x sn\k x + A y sink y ). (A18) 

Of course, the order parameters, P, Q, IT, and A, should be 
determined self-consistently by satisfying the following con- 
ditions: 



Px — (^j q 7i+:e,ck) 5 Qx — (^ia'n+:r,a) j 

ITrr = (^n a ^3,i— £,a) — {h\i a hi,i- x ,a) 

= (^2ia^34+i,a) = (^ia^i+^a) > 
A x = (/lliJ./l3,i-«,T) — (/lH1^4,l-*,t) 
= (/l2iJ./l3,i+S,T) = (^2il^4,i+«,T)) 



(A19) 



and similar equations for the y-direction. 

As seen in Eq. (IA14b and subsequent equations, Ck, dk, ek, 
and /k are in general complex. In this case it is not easy to 
obtain analytic expressions for the ft,-fermion spectrum. For- 
tunately, however, it turns out that we can make the following 
approximations: 



Ck — cm, dk = dm, 
eik — e_Rk, e-2k — em, 
fik — fm, /2k — fsk- 



(A20) 



The justification for the above approximation will be given in 
Appendix I A 3 1 where we discuss the physical meaning of the 
approximation and provide supporting numerical results. For 
the time being, however, we proceed to derive saddle-point 
equations by accepting Eq. ( 1A201 i. 



2. Self-consistent saddle-point equations 

After the Bogoliubov transformation, the saddle-point 
Hamiltonian in Eq. dA8t can be written in the following way: 



H = Ne + ^2J2[b k -Qi(k)} 

k 1 = 1 

4 

+ £]>>(k) 7 t Q (k) 7;a (k) 
k 1=1 

+ E { w vi( k )vL(k)Vi*(k) + ^Ki(k) - C k ] J 

+ \ w MvL(^)m a (k) + f K 2 (k) - c k ] 



(A21) 



where 



w r , 1 (k) = ^f(C k + A k ) 2 -B 2 , 
w v2 (k) = J(Ck- A k )*-B*, 



fii(k) = V («k - b k + c k - dk) 2 , 
n 2 (k) = V ( a k -b k -c k + 4) 2 , 



JMk) = y/(a k + b k -c k - d k ) 2 + |e k + / k | 2 , 

4 (k) = V / («k + 6k + c k + d k ) 2 + |e k + /k| 2 , (A22) 

and we set c k = c Rk , d k = d Bk , ek = e_R k , and / k = 
fm- I n the above equations, 7/ lQ ,(k) and r] 2a {k) describe 
two bosonic particles derived from the linear combination of 
t a (k) and r a (k). Similarly, 7;(k) denote four fermionic Bo- 
goliubov quasi-particles obtained from h m (k) and h] n (— k) 
[to- 1,2,3,4]. 

Now the ground state energy per unit cell is given by 



(g)gr 

N 



k 1=1 

+ jf E \ Kl(k) _ Ck] + ^ E I K^ k ) - <K ■ 

k ~ 



N ^ 2 

k 



(A23) 



The four saddle-point equations for s, d, £, and /i can be ob- 
tained by minimizing the ground state energy as follows: 



de 



d(e gr /J) de gr d(e gr /J) 



dd 2 







(A24) 



where 



de 



9 T _ c 



= 5-2F 



2d 2 - Ax 



3 J_ 

27V 



E 



C k + A k , Ck-A k 



W,,i(k) Wr,2(k) 



8 



d(e gr /J) _ 3 j 
9s 2 2 J 



iV 



^(cos k x + cos fcy) 



Qk + h. - Ck - d k 

r» 3 (k) 

«k + h + c k + 4 



\/2 ^ — tt- (c°s + cos fcj, ) 



n 4 (k) 

ek + /k e k + f k 

r> 3 (k) n 4 (k) 



<9e 



-A{x -<?-!) 



«k - bk - c k + 4. a k - fe k + c k - c? k 



fii(k) 



n 2 (k) 



ak + b k - c k - d k a k + b k + c k + d k 



n 3 (k) 



Q 4 (k) 



d(e gr /J) , A V, 9 V 

dd? ~ j j j 



2^(cosfc x + cos fc y ) 

ak + b k + c k + d k 



N 

2t'/J 
N 



J^(cos k x + cos k y ) 



+ V^i^— ^(COS^ + COS&„) 



fi 4 (k) 

Ck + /k . e k + / k 



n 3 (k) n 4 (k) 

(A25) 



E(k)/f 




(0,0)(^/2,7i/2)(0,7t) (0,0) 

(kx.ky) 

FIG. 6: (Color online) Band structure of the ft-fermion at half- 
filling, (a) The full band structure of the mean-field Hamiltonian, 



= H, 



(") 



+ #4 



(b) The band structure 



obtained if Hfermion is approximated to be H t 



(0) 



#4 



(C) 



The additional band structure due to the remaining term, 



For convenience, we plot the band structure of H^ nion + , |_ 



terms of the bonding and anti-bonding fermionic operators 
that can be written as follows: 



Vi+a 



1 

71 
i 

71 
i 

71 



(h 



Ha — Il2ia) 



3ia '»£ia) 



(A27) 



Self-consistency conditions for eight order parameters can 
also be written in terms of new quasi-particle spectra: 



P, 

Q, 
n, 



3 1 

In 



cos k x 



C k + A k C k - A k 



3 1 

an 



cos k x 



W v l(k) Wr,2(k) 

Bk B k 



--^cosfc, 



w 71 i(k) w v2 (k) 
a-k + bk - c k - d k 



AN 

1 1 

47V 



cos k a 



n 3 (k) 

a k + b k + c k + 4 



A. 



1 1 



cos k x 



ft 4 (k) 

ek + /k e k + /k 
fi 3 (k) rj 4 (k) 



(A26) 



and additional four self-consistency conditions for P y , Q y , 
Hy, and A y are obtained when cos k x is replaced by cos k y . 



3. Validity of the approximated /i-fermion dispersion 



In order to justify the approximations used in Eq. (|A20| i, we 
start by expressing the fermionic part of the Hamiltonian in 



where h+ and /i_ describe the bonding and anti-bonding state 
for the horizontal dimers, respectively. v+ and t>_ play exactly 
the same role for the vertical dimers. 

The fermionic part of the Hamiltonian can be decomposed 
into three parts: 



TT _ rr(0) , TT 

-nfermion — -"fermion 4 



(A28) 



where 



(o) 

fermion 



= 5^ 0k [4 a (k)ft +a (k)-ftl (k)ft_ (k) 

k 

+ 4 (kK„(k)-«I a (k)i;_ (k) 
J2bk[hl a (k)h +a (k) + hl a (k)h_ a (k) 

k 

+ vUk)v +a (k) + vl a (k)v_ a (k) 



(A29) 
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H 



++, — 



= CRk [ h +a(^>+a(k) + vl a {k)h +a (k) 
k 

+ hl a (k)v^ a (k)+vl a (k)h^ a (k) 

+ d ™ [ h +a(k)v+a(k) + vl a {k)h +a (k) 
k 

-ftl a (k)t;_a(k)-t;i a (k)/i_«(k) 
+ E {eflk[«t T (k)^i(-k) + ^ T (k)< x (-k) 

k 

+ vl T (k)/il ; (-k) + ^ T (k)u^(-k)] +H. c.} 

+ E {.m^t^V-*) + ^l T ( k Ki(- k ) 

k 

- ^ T (k)/iLx(-k) - ^(kjvl^-k)] +H. c.}, 

(A30) 

and 

#+-,-+ 

= J2^[ V -a( k ) h +a(k) - hl a (k)v„ a (k) 
k 

+vi a (k)h_ a (k)-hi a (k) v+a (k) 

+ J2idi k [vUk)h_ a (k) - hl a (k)v +a (k) 

k 

+ hl a (k)v^ a (k)-vl a (k)h +a (k) 
+ E {ie Jk [^ T (k)/i^(-k) - ^ T (k>l ; (-k) 

k 

+ ^ T (k)/il ; (-k) - ftl t (k)vj. i (-k)] +H. c.} 

+ E {*//k[4 T (k)^ x (-k) - /ii T (k)< x (-k) 

k 

+ ^ T (k)^(-k) - wi T (k)^(-k)] + H. c.}, 

(A31) 



where it is important to note that describes the 
fermion hopping and pairing processes within the band of the 
same bonding character. On the other hand, t |- con- 
tains mixing processes between the bonding and anti-bonding 
bands. The approximation used in Eq. dA20b corresponds to 
the omission of i |_ from H f er mion- 



To judge the validity of this approximation, let us first con- 
sider the effect of hopping terms alone from the Hamiltonian. 
In Fig. [6] (a) we plot the full band structure of -fff er mion at half 
filling. Since the unit cell consists of four spins, there are 
four bands. The two bands lying high in energy are almost 
flat while the other two are dispersive. In Fig. 0(b) we show 
the band structure obtained when Hf am { oa is approximted to 



be H, 



(o) 

term; 



H ++ , As one can see, the difference be- 



tween the full and the approximated band structure is minor. 
The almost flat two bands in the full Hamiltonian, which are 
high in energy and so do not play important roles at small 
hole doping, become completely localized. The shapes of the 
other two dispersive bands are only slightly modified. Now, 
to explicitly demonstrate the minor influence of the remaining 

part, i.e., . |_, we consider an extreme situation where 

the role of i |_ is maximized. That is, we consider the 

band structure of H^ mion + . |_. As seen in Fig. |6](c), 

the resulting band structure is almost featureless. Thus, the 
main effect of i |_ would be providing only a weak dis- 
tortion to the bands of H ++: 



Since normal mixing between the bands with different 
bonding characters generates minor effects, it is reasonable 
to assume that the pairing process between different bonding 
bands also does not modify the band structure significantly. 
It can be further argued that these results are valid also for 
the case of finite doping because the fermionic band structure 
has similar structure to the half-fillied case as can be seen in 
Eq. (lA22b and there is no level crossing upon doping. In con- 
clusion, we believe that the approximation used in Eq.( lA201 > 
is accurate enough to capture the essential physics while it 
renders convenient simplification in calculations. 
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